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A pulse tube type Stirling machine with warm gas-driven displacer which has a displacer rod is discussed 
with numerical simulation when it is used as a cryogenic refrigerator, room temperature refrigerator and 
engine. It has both the advantages of gas-driven-stirling machine with high efficiency and simplicity and 
the advantages of pulse tube machine with no moving parts at low temperatures. A nodal analysis 
method that includes the linear motor and the displacer in the machine is introduced. Numerical results 
show that it has high potential to be used as the cryogenic refrigerator, room temperature refrigerator 
and engine. In this type of machine, there is an optimum phase angle between displacer and piston, 
and an optimum swept volume ratio of displacer over compressor for efficiency. The phase angle and 
swept volume ratio can be adjusted by the natural frequency of the displacer and the diameter of the dis¬ 
placer rod when it is used as a refrigerator. 

© 2010 Elsevier Ltd. All rights reserved. 


1. Introduction 

The advantage of the common pulse tube refrigerators, such as 
orifice, double inlet, and inertance tube pulse tube refrigerators, is 
that there is no moving part at low temperatures and only one 
moving part at room temperature [1-5]. The expansion work in 
the pulse tube of the refrigerator is changed to heat, therefore 
the theoretical efficiency is lower compared to a Carnot refrigera¬ 
tor. In a warm piston or warm displacer type pulse tube refrigera¬ 
tor [6-11] besides a compressor as moving part at room 
temperature, there is a piston or a displacer at the warm end of 
the pulse tube to recover the expansion work. Though it is com¬ 
plex, the theoretical efficiency is the same as that of a Carnot refrig¬ 
erator if all else is perfect, and any phase angle between the 
pressure and the flow can be achieved to minimize the regenerator 
loss [8]. In this paper, a nodal analysis numerical method was 
introduced for a pulse tube machine with a warm gas-driven dis¬ 
placer which has a displacer rod [11], and the machine was simu¬ 
lated by a numerical method when it is used as cryogenic 
refrigerator, room temperature refrigerator, and engine. The 
numerical method is similar to that in Ref. [5], and the motions 
of the linear motor and the displacer are included. 


2. Structure of the novel machine 

Fig. 1 shows the schematic of the pulse tube type Stirling ma¬ 
chine that has a displacer with a displacer rod as a phase shifter. 
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It includes a compressor, a displacer part and a fixed part. The fixed 
part includes a cooler, a regenerator, an absorption heat exchanger, 
a pulse tube, a connecting tube, and gas supply tube. The compres¬ 
sor includes a motor house, a linear motor, a group of flexible bear¬ 
ings and a piston. The displacer part includes a buffer, a group of 
flexible bearings, a displacer and a displacer rod. 

The displacer with the rod is suspended by flexible bearing to 
form a clearance seal. The displacer with the rod and the flexible 
bearings form an oscillation system. The displacer oscillates due 
to the pressure oscillation in working spaces and the buffer. The 
expansion power of the pulse tube is obtained on the displacer 
expansion space and is transferred to the displacer back space that 
joins the compression space to compress the working gas. 

When it is used as a refrigerator, cooling power is obtained from 
the absorption heat exchanger which is at a low temperature. The 
electric power is input to the linear motor. Compared to the iner¬ 
tance tube pulse tube refrigerator [5], this machine has a higher 
efficiency because the expansion power from the pulse tube is 
recovered by the displacer. When operating as a room temperature 
refrigerator, the advantage is even larger because the expansion 
power from the pulse tube is not so small compared to the input 
power of the compressor. 

When it is used as an engine, the heat is supplied by the absorp¬ 
tion heat exchanger which is kept at a high temperature. The linear 
motor becomes a linear generator from which electric power can 
be obtained. Compared to a Stirling engine, the advantage of this 
type of machine is that there is no moving part at high tempera¬ 
ture, and the moving parts are at room temperature and can be 
separated from the fixed parts. 

In Fig. 1, the displacer is driven by gas pressure difference. Only 
the linear motor needs a sinusoidal wave current which can be 
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Nomenclature 



a 

swept volume ratio 

R 

gas constant 

A 

cross gas flow area 

Pm 

motor coil resistance 

As 

cross solid area 

S 

enthalpy 

Al 

heat transfer area per meter 

S 

enthalpy flow 

B 

motor force factor 

t 

time 

C v 

constant volume specific heat of gas 

T 

temperature of gas 

Cp 

constant pressure specific heat of gas 

To 

room temperature 

Cs 

specific heat of solid 

T s 

temperature of matrix 

D e 

hydro-diameter 

u 

speed 

D P 

diameter of piston 

U 

voltage 

Dd 

diameter of displacer 

U 0 

voltage amplitude 

D r 

diameter of displacer rod 

Vc 

volume of compressor 

E 

exergy flow 

Vcs 

dead volume of compressor 

fr 

friction coefficient 

Vq 

swept volume of compressor 

fp 

resistance of piston 

V e 

displacer volume 

fd 

resistance of displacer 

Ves 

displacer dead volume 

Fp 

piston clearance leakage factor 

Ve 0 

displacer swept volume 

F d 

rod clearance seal leakage factor 

VeB 

displacer back volume 

H 

enthalpy flow 

w 

compressor power 

Hpj 

pulse tube enthalpy flow 

w m 

motor input power 

HrEG 

regenerator loss 

X 

coordinator along the length 

i 

current 

Xp 

displacement of piston 

I< P 

stiffness of piston spring 

x d 

displacement of displacer 

K d 

stiffness of displacer spring 

X Po 

distance between piston head and cylinder head at 

L 

inductance of linear motor 


standstill 

m 

mass flow rate 

*d 0 

distance between displacer head and displacer cylinder 

ifip 

gas leakage from piston seal 


head at standstill 

m d 

gas leakage from rod seal 



m s 

mass of solid 

Greek letters 

m p 

mass of piston and mover of linear motor 

a 

heat transfer coefficient 

m d 

mass of displacer 

rj 

engine efficiency 

P 

pressure 

P 

density of gas 

Pp 

pressure in compressor cylinder 

Ps 

density of matrix 

Pd 

pressure in displacer cylinder 

T 

period 

Pb p 

pressure in motor house 


phase angle between compressor and displacer 

P Bd 

pressure in buffer 

CO 

angle frequency 

a 

cooling power or heating power 




provided by an AC power source or a simple inverter. There is no 
moving part at low temperature or high temperature range. In this 
sense the machine has the same advantage of both a pulse tube 
machine and a gas driven-stirling machine. We can call this ma- 



Fig. 1 . Schematic drawing of the machine: 11. cooler; 12. regenerator; 13. 
adsorption heat exchanger; 14. pulse tube; 15. connecting tube; 16. gas supply 
tube; 21. compressor spring; 22. linear motor; 23. motor house; 24. piston; 25. 
compression space; 31. displacer spring; 32. rod; 33. buffer; 34. displacer; 35a. 
expansion space; 35b. back space. 


chine a “pulse tube Stirling machine with warm gas-driven displac¬ 
er” or, alternatively, a “Stirling machine with pulse tube”. 

In Fig. 1, the buffer does not consume any work theoretically if 
the gas spring loss, friction loss, and gas leakage loss from the 
clearance seal around displacer rod are not considered. 

3. Numerical simulation method 

The numerical simulation is separated into two steps. In step 
one, the numerical model does not include the motion of the linear 
motor and the displacer. The motions of the displacer and the pis¬ 
ton are assumed as sinusoidal to investigate the basic character. In 
step two, the numerical model includes the motion of the linear 
motor and the displacer to investigate how the displacer moves 
in the cryogenic refrigerator. 

Assuming that the volume of the gas supply tube connecting the 
back space of the displacer and compression space of the compres¬ 
sor and the cooler is neglected, the compression space and the back 
space of the displacer are treated as single volume. 

3.1. Sinusoidal displacement of piston and displacer 

The numerical simulation was done by a nodal analysis method 
based on Ref. [5], and was mainly based on the assumptions that 
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the system is one dimensional, and that the working gas is an ideal 
gas. The governing equations are as follows. 

Continuity equation 


d d . 

-(pA)+-(puA) 


0 


Momentum equation 


dP f r 1 2 u 

dx + D e 2 pU |u| 


0 


Energy equation for gas 

l, P AC,Tnl^AC,T ) + aA I ,T-T, l+ ^ 

Energy equation for matrix 

dT s 


0 


AsP s C. 


dt 


ocA l (T-T s ) 


Ideal gas state equation 
P = pRT 

Compression volume of the compressor 
Vc = V cs + 0.5 V Co (1 + cos (cot - q>)) 

Expansion volume of the displacer 
V E = Ves + 0.5V £o (1 + cos(cot)) 

Back volume of the displacer 
Veb — Ves + 0.5V£ 0 (1 — cos (cot)) 

Swept volume ratio of the displacer over the compressor 


a = 


V 


E o 


Vc n 


Compression power of the compressor 


W 


U 


PdV 


c 


T J 

Enthalpy flow 

H = i j) mC p Tdt 

Entropy flow 

5 = 1 £ msdt 

Exergy flow 
E = H -ST o 
Cooling power 

Q. = ffpT — Hreg 

COP of refrigerator 
COP = Q/W 

Efficiency of engine 
rj = W/Q 


( 1 ) 


( 2 ) 


( 3 ) 


( 4 ) 


( 5 ) 


( 6 ) 


( 7 ) 


( 8 ) 


( 9 ) 


( 10 ) 


( 11 ) 


( 12 ) 


(13) 


(14) 


(15) 


(16) 


The equivalent P-x diagram in the pulse tube definition is the same 
as that in Ref. [5]. 

3.2. Equations of motion for piston and displacer 

As shown in Fig. 1, the linear motor drives the piston oscillation. 
The electric force of the linear motor is from the electric current. 
The forces on the piston are from the motor, the flexible spring, 


the friction and the gas pressure difference on the two ends of 
the piston. The forces on the displacer are from the flexible spring, 
friction, the gas pressure difference between the buffer and the 
pulse tube on the rod, and the gas pressure difference on two ends 
of the displacer. Assuming that the friction factor of the clearance 
is constant for the gas leakage through the clearance seal, and that 
the electric force of the linear motor is proportional to the electric 
current, the following equations are obtained. 

Differential equation of the motion of the piston 

iripXp T fpXp T kpXp — Bi -(- — Dp (Pgp Pp ) (17) 

Differential equation of the motion of the displacer 

m d x d +f d x d + k d x d = ^(D 2 d - D 2 )(P p - P d ) +^D 2 r (P Bd - P d ) (18) 

Equation of electrodynamic driver 
U = Uo sin (cot) = Bx p + iR m +L^l (19) 


The gas leakage through the clearance seal between the com¬ 
pression volume and the motor house 



The gas leakage through the clearance seal between the back 
volume and the buffer 



Motor input power 





Uidt 




With the above equations, the simulation of the displacement of the 
piston and the displacer can be achieved. Then the volume V c , V E , 
and V EB in Eqs. (6)—(8) can be replaced by the following equations. 
Compression volume of the compressor 


= ! d p( x p„ -Xp) 

(23) 

Expansion volume of the displacer 


= ^D 2 d (x do -x d ) 

(24) 

Back volume of the displacer 


: =E (D 2 -D 2 )(x do +x d ) 

(25) 


Table 1 

Calculation condition of cryogenic refrigerator. 


Regenerator 

0120 mm x 50 mm, porosity 0.7, mesh 
wire diameter 0.04 mm 

Pulse tube 

040 mm x 500 mm 

Cooler height, width, length 

0.33 mm x 2837.25 mm x 40 mm 

Absorption heat exchanger 
height, width, length 

0.5 mm x 829.38 mm x 40 mm 

Compressor swept volume 

500 cm 3 

Refrigeration temperature 

77 K 

Room temperature 

300 K 

Frequency 

50 Hz 

Charge pressure of helium gas 

2.5 MPa 
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Fig. 2. (a) COP vs. phase angle and swept volume ratio, (b) Cooling power vs. phase angle and swept volume ratio, (c) Input power vs. phase angle and swept volume ratio, (d) 
Enthalpy flow of pulse tube vs. phase angle and swept volume ratio, (e) Enthalpy flow of regenerator vs. phase angle and swept volume ratio, (f) Pressure ration in pulse tube 
vs. phase angle and swept volume ratio. 
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4. Numerical results with sinusoidal displacement 

In this section, the numerical simulation with sinusoidal dis¬ 
placement of compressor and displacer is conducted when the 
pulse tube Stirling machine is used as cryogenic refrigerator, as 
room temperature refrigerator, or as engine. 

4A. Cryogenic refrigerator 

In this section, a kilowatt order refrigerator was simulated. 
With given fixed parts, the swept volume ratio of the displacer to 
the compressor, the phase angle between the compressor and the 
displacer are investigated. This is the basic information required 
for the design of this type of refrigerator. 

Table 1 shows the basic data of the refrigerator for the simula¬ 
tion. With this condition, Fig. 2 shows the COP, cooling power, in¬ 
put power, pressure ratio, enthalpy flow of the pulse tube and 
enthalpy flow of the regenerator change with the swept volume ra¬ 
tio of the displacer to the compressor and the phase angle between 
the compressor and the displacer. In Fig. 2, the swept volume ratios 
are 0.2, 0.4, 0.6, 0.8 and 1. There is an optimum phase angle for 
each of COP, cooling power, input power, and the enthalpy flow 




of the pulse tube. The optimum phase angle for cooling power is 
larger than that for the COP, and is smaller than that for the input 
power. The enthalpy flow of the regenerator increases with the in¬ 
crease of the phase angle. There is an optimum swept volume ratio 
of the displacer to the compressor for COP and cooling power. In 
the calculation range, the input power and the enthalpy flow of 
the pulse tube increase with the increase of the volume ratio of 
the displacer to the compressor. When the phase angle is smaller 
than 50°, there is an optimum swept volume ratio for minimum 
enthalpy flow of the refrigerator. When the phase angle is over 
50°, the enthalpy flow of the regenerator increases with the in¬ 
crease of the swept volume ratio of the displacer to the compres¬ 
sor. The pressure ratio decreases with the increase of the phase 
angle and increases with the increase of the volume ratio of the 
displacer to the compressor. 

It should be mentioned that the COP is over 0.2, which indicates 
this machine could get higher efficiency than that of the inertance 
tube pulse tube refrigerator [5]. One of the reasons is that the 
expansion power of the pulse tube is recovered by the displacer. 
Another reason may be the high pressure ratio. At optimum COP 
point, the pressure ratio is 1.5. High pressure ratio also means this 
type of refrigerator has high power density. 




Fig. 3. (a) P-x diagrams in pulse tube, (b) Mass flow rate and pressure wave at the cold end of regenerator, (c) PV diagrams of displacer and compressor, (d) Enthalpy flow, 
exergy flow and temperature distribution in the machine. 
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Fig. 3 shows the P-x diagrams in the pulse tube, the pressure 
wave and the mass flow rate at the cold end of the pulse tube, 
the PV diagrams of compressor and displacer, enthalpy flow, exer- 
gy flow, and temperature distribution of the machine at optimum 
COP point when the phase angle is 20° and swept volume ratio is 
0.6. There is a long distance between the P-x diagrams, which indi¬ 
cates that the pulse tube length is large enough. Sufficient distance 
between the P-x diagrams at the cold end and warn end of the 
pulse tube means there is a sufficient separation between cold 
gas and warm gas. If this were not the case, cold gas and warm 
gas flowing to the same place may generate large shuttle loss in 
the pulse tube. The pressure is ahead of a large phase angle to mass 
flow rate, which indicates high efficiency condition of the regener¬ 
ator. The PV diagram of the displacer is not so small compared to 
the PV diagram of the compressor, which indicates that the gain 
from the displacer is not so small. In Fig. 3d, the flat line of enthal¬ 



py flow indicates the regenerator and the pulse tube. The flat line 
of temperature indicates the cooler, absorption heat exchanger, 
and connecting tube. The exergy flow decreases from left to right 
and remains almost constant in the pulse tube to the hot end 
where it is taken off by the displacer. 

Fig. 4 shows the COP and cooling power changes with swept 
volume of the compressor under the condition of 20° of the phase 
angle between the compressor and the displacer and 0.6 of swept 
volume ratio of the displacer over the compressor. In the calcula¬ 
tion range, with the increase of the compressor swept volume, 
the cooling power, input power, pressure ratio increase, and there 
is an optimum point for COP. 



Refrigeration temperature, K 

Fig. 5. Cooling power, input power and COP of room temperature refrigerator. 


Table 2 

Calculation condition of room temperature refrigerator. 



Regenerator 
Pulse tube 

Cooler height, width, length 
Absorption heat exchanger 
height, width, length 
Compressor swept volume 
Displacer swept volume 
Phase angle 
Room temperature 
Frequency 

Charge pressure of helium gas 


030 mm x 30 mm, porosity 0.7449, mesh 

wire diameter 0.05 mm 

020 mm x 120 mm 

0.33 mm x 353.43 mm x 20 mm 

0.33 mm x 353.43 mm x 20 mm 

15 cm 3 
15 cm 3 
30° 

333.15 K 
50 Hz 
3 MPa 


Table 3 

Calculation condition of engine. 


Fig. 4. (a) COP and cooling power vs. swept volume of compressor, (b) Input power 
and pressure ratio vs. swept volume of compressor. 


Regenerator 
Pulse tube 

Cooler height, width, length 
Absorption heat exchanger 
height, width, length 
Compressor swept volume 
Heater temperature 
Room temperature 
Frequency 

Charge pressure of helium gas 


030 mm x 60 mm, porosity 0.7, mesh 

wire diameter 0.1 mm 

030 mm x 200 mm 

0.33 mm x 353.43 mm x 40 mm 

0.33 mm x 353.43 mm x 40 mm 

50 cm 3 
900 K 
300 K 
100 Hz 
3 MPa 
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Fig. 6. (a) Efficiency vs. swept volume ratio and phase angle, (b) Input heat vs. swept volume ratio and phase angle, (c) Output power vs. swept volume ratio and phase angle, 
(d) Enthalpy flow of pulse tube vs. swept volume ratio and phase angle, (e) Enthalpy flow of regenerator vs. swept volume ratio and phase angle, (f) Pressure ratio in pulse 
tube vs. swept volume ratio and phase angle. 
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The above results show that there is an optimum swept volume 
ratio of the displacer to the compressor and an optimum phase an¬ 
gle between the compressor and the displacer for COP at a given 
swept volume of the compressor. 

4.2. Room temperature refrigerator 

Fig. 5 shows the input power, cooling power and COP change 
with the refrigeration temperature under the conditions shown 
in Table 2. The COP and cooling power increase with the increase 
in refrigeration temperature, and the input power decreases with 
the increase in refrigeration temperature. Near room temperature 
the COP is larger than 2. A COP of 2 can not be achieved by the iner- 
tance tube pulse tube refrigerator because its theoretical efficiency 
is less than unity. 

In Fig. 5, the COP is 7 when the temperature of the absorption 
heat exchanger is the same as the room temperature. Theoretically, 
the machine works as an engine when the temperature of the 
absorption heat exchanger is higher than that of the cooler. The 
machine generates an output power of 7 W when the temperature 
of the absorption heat exchanger is 393.2 K. From the room tem¬ 




perature to a point where the temperature of the absorption heat 
exchanger is higher than the room temperature and the output 
power is zero, the machine just consumes the input power, or 
the machine cannot move due to loss. 

When the temperature of the absorption heat exchanger is the 
same as the room temperature and the temperature of the cooler is 
higher than the room temperature, the machine works as a heat 
pump. 

This type of machine can be considered as the merge of pulse 
tube and Stirling refrigerator. So it can be very small and be oper¬ 
ated at a wide range of refrigeration temperatures. 

4.3. Engine 

In this section, a kilowatt order engine was simulated. With gi¬ 
ven fixed parts, the swept volume ratio of the displacer to the com¬ 
pressor, the phase angle between the compressor and the displacer 
are investigated. The calculation condition is shown in Table 3. 
Flere, the higher frequency is used for higher power density. 

Fig. 6 shows the input heat, efficiency, output power, pressure 
ratio, enthalpy flow of the pulse tube, and enthalpy flow of the 




Fig. 7. (a) P-x diagrams in pulse tube, (b) Mass flow rate and pressure wave at the hot end of the pulse tube, (c) PV diagrams of displacer and compressor.(d) Enthalpy flow, 
exergy flow and temperature distribution in the machine. 
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Table 4 

Calculation condition of cryogenic refrigerator. 


Regenerator 

0120 mm x 50 mm, porosity 0.7, mesh 
wire diameter 0.04 mm 

Pulse tube 

040 mm x 500 mm 

Cooler height, width, length 

0.33 mm x 2837.25 mm x 40 mm 

Absorption heat exchanger 
height, width, length 

0.5 mm x 829.38 mm x 40 mm 

Refrigeration temperature 

77 K 

Room temperature 

300 K 

Frequency 

50 Hz 

Charge pressure of helium gas 

2.5 MPa 

Voltage 

283 V 

Piston diameter 

150 mm 

Piston spring stiffness 

500 x 10 3 N/m 

Piston weight 

11.5 kg 

Displacer diameter 

120 mm 

Displacer spring stiffness 

50 x 10 3 N/m 

Displacer weight 

0.5 kg 

Displacer rod diameter 

30 mm 



regenerator change with the swept volume ratio of the displacer to 
the compressor and the phase angle between the compressor and 
the displacer. In the figure, the swept volume ratios are 0.2, 0.4, 
0.6, 0.8 and 1. There is an optimum phase angle for efficiency, input 
heat, output power, and the enthalpy flow of the pulse tube. The 
optimum phase angle for output power is larger than that for effi¬ 
ciency, and is smaller than that of input heat. The enthalpy flow of 
the regenerator increases with the increase of the phase angle. 
There is an optimum swept volume ratio of the displacer over 
the compressor for efficiency. In the calculation range, the input 
heat, the output power and the enthalpy flow rate of the pulse tube 
increase with the increase of the volume ratio of the displacer to 
the compressor. When the phase angle is smaller than 70°, there 
is an optimum swept volume ratio for minimum enthalpy flow of 
the refrigerator. When the phase angle is over 70°, the enthalpy 
flow of the regenerator increases with the increase of the swept 
volume ratio of the displacer to the compressor. The pressure ratio 




Fig. 8. (a) Voltage and current, (b) Displacement of displacer and piston. (cO Pressure in compressor, displacer, motor house and displacer buffer. 
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Fig. 9. (a) COP vs. piston weight, displacer weight and rod diameter: 1. rod 30 mm, piston 11 kg; 2. rod 30 mm, piston 11.5 kg; 3. rod 30 mm, piston 12 kg; 4. rod 40 mm, 
piston 11 kg; 5. rod 40 mm, piston 11.5 kg; 6. rod 40 mm, piston 12 kg; 7. rod 50 mm, piston 11 kg; 8. rod 50 mm, piston 11.5 kg; 9. rod 50 mm, piston 12 kg. (b) Cooling 
power vs. piston weight, displacer weight and rod diameter: 1. rod 30 mm, piston 11 kg; 2. rod 30 mm, piston 11.5 kg; 3. rod 30 mm, piston 12 kg; 4. rod 40 mm, piston 11 kg; 
5. rod 40 mm, piston 11.5 kg; 6. rod 40 mm, piston 12 kg; 7. rod 50 mm, piston 11 kg; 8. rod 50 mm, piston 11.5 kg; 9. rod 50 mm, piston 12 kg. (c) Amplitude of displacement 
of displacer vs. piston weight, displacer weight and rod diameter: 1. rod 30 mm, piston 11 kg; 2. rod 30 mm, piston 11.5 kg; 3. rod 30 mm, piston 12 kg; 4. rod 40 mm, piston 
11 kg; 5. rod 40 mm, piston 11.5 kg; 6. rod 40 mm, piston 12 kg; 7. rod 50 mm, piston 11 kg; 8. rod 50 mm, piston 11.5 kg; 9. rod 50 mm, piston 12 kg. (d) Ratio of 
displacement vs. piston weight, displacer weight and rod diameter: 1. rod 30 mm, piston 11 kg; 2. rod 30 mm, piston 11.5 kg; 3. rod 30 mm, piston 12 kg; 4. rod 40 mm, piston 
11 kg; 5. rod 40 mm, piston 11.5 kg; 6. rod 40 mm, piston 12 kg; 7. rod 50 mm, piston 11 kg; 8. rod 50 mm, piston 11.5 kg; 9. rod 50 mm, piston 12 kg. (e) Phase angle between 
piston and displacer vs. piston weight, displacer weight and rod diameter: 1. rod 30 mm, piston 11 kg; 2. rod 30 mm, piston 11.5 kg; 3. rod 30 mm, piston 12 kg; 4. rod 40 mm, 
piston 11 kg; 5. rod 40 mm, piston 11.5 kg; 6. rod 40 mm, piston 12 kg; 7. rod 50 mm, piston 11 kg; 8. rod 50 mm, piston 11.5 kg; 9. rod 50 mm, piston 12 kg. 
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increases with the phase angle increase, and decreases with the in¬ 
crease in the swept volume ratio of displacer to compressor, which 
is opposite to that in the refrigerator. 

Fig. 7 shows the P-x diagrams in the pulse tube, the pressure 
wave and mass flow rate at the hot end of the pulse tube, the PV 
diagrams of compressor and displacer, enthalpy flow, exergy flow, 
and temperature distribution in the machine at swept volume ratio 
of 0.6 and phase angle of 20°. The pressure is delayed by a large 
phase angle with respect to the mass flow rate. There is a sufficient 
distance between the P-x diagrams. The PV diagram of displacer is 
larger than that of compressor. In Fig. 7d, the flat line of enthalpy 
flow indicates the regenerator and the pulse tube. The flat line of 
temperature indicates the cooler, absorption heat exchanger, and 
connecting tube. The enthalpy loss in the regenerator is from the 
hot end to the warm end, so it is negative. The exergy flow in the 
heater increases, and keeps nearly constant in the pulse tube. 

The above results show that there is an optimum phase angle 
and swept volume ratio for the efficiency of the engine at a fixed 
compressor swept volume. This is the same when it is used as 
the refrigerator. 

5. Numerical results with linear motor and displacer 

In the following section, the mass of the displacer, the diameter 
of the rod, and the mass of the piston are discussed when the ma¬ 
chine is used as the cryogenic refrigerator under the condition that 
the linear motor and displacer are included in the numerical mod¬ 
el. Table 4 shows the calculation condition. 

Fig. 8 shows the voltage, current, displacement of piston and 
displacer, and pressure with the condition shown in Table 4. The 
amplitudes of the displacement of the piston and the displacer 
are 16 mm and 13 mm, respectively. The phase angle between 
the displacement of the displacer and the piston is 24°. The phase 
angle between the current and the displacement of the piston is 
86°, which means that the current and the velocity of the piston 
are almost in phase, and the motor oscillates near its resonance 
point. The pressure in the motor house and the buffer is near the 
average level. 

Fig. 9 shows the COP, cooling power, amplitude of displacement 
of the displacer, amplitude ratio of the displacer to the piston, 
phase angle between the displacer and the piston change with 
the mass of the displacer, the diameter of the displacer rod, and 
the mass of the piston. For each displacer rod diameter, there is 
an optimum displacer mass for the COP, and the optimum point 
is hardly influenced by the piston mass. The cooling power in¬ 
creases with the decrease of the displacer mass and with the in¬ 
crease of the piston mass. The amplitude of the displacement of 
the displacer increases with the increase of the displacer mass 
and the piston mass. The amplitude ratio of the displacer to the 
piston has little relation with the mass of the piston, and increases 
nearly in a linear fashion with the increase of the mass of the dis¬ 
placer. The phase angle between the displacer and the piston is 
scarcely influenced by the mass of the piston, and almost decreases 
linearly with the decrease of the mass of the displacer. 

At optimum COP point, for rods of 30 mm, 40 mm, and 50 mm, 
the masses of the displacer are 0.5 kg, 0.75 kg, and 1.1 kg, respec¬ 
tively, the phase angle between the displacer and the piston is 


24°, 24°, and 22°, respectively, and the ratios of the amplitude of 
the displacement of the displacer to the piston are 0.8, 0.86, and 
0.95, respectively. This means that a bigger rod diameter gives a 
bigger amplitude ratio of the displacer to the piston at the opti¬ 
mum point. So the bigger rod diameter has a higher swept volume 
ratio of the displacer to the piston. 

In Fig. 9, the increase of the mass of the displacer means that the 
natural frequency of the displacer is decreased. So the natural fre¬ 
quency is a strong parameter to adjust the phase angle between 
the piston and the displacer. 

In this section, the numerical results show that the phase angle 
between displacer and piston can be controlled by the mass of the 
displacer, and the displacement ratio of the displacer to the piston 
can be controlled by the rod diameter and the mass of the 
displacer. 

6. Conclusion 

We have introduced a nodal analysis numerical model for the 
simulation of the pulse tube Stirling machine with warm gas-dri¬ 
ven displacer and a displacer rod. The model was use to simulate 
a machine when it is used as cryogenic refrigerator, room temper¬ 
ature refrigerator and engine. The results are useful for under¬ 
standing this type of machine. For given fixed parts and a given 
swept volume, there is an optimum phase angle between the com¬ 
pressor and displacer for COP when it is used as refrigerator, and 
for efficiency when it is used as engine. When it is used as cryo¬ 
genic refrigerator, the mass of the displacer can be used to adjust 
the phase angle between the compressor and the displacer, and 
the amplitude ratio of the displacer to the piston when the rod 
diameter is fixed. At the optimum COP point, the larger rod diam¬ 
eter has the higher amplitude ratio of the displacer to the piston. 
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